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A JET AND A FASTER MOVING COAXIAL STREAM 
by Leo F. Donovan 
Lewis Research Center 

SUMMARY 

An analysis of turbulent mixing of a jet with a faster moving coaxial stream is pre- 
sented, and a similarity solution is obtained for large initial ratios of coaxial- stream ve- 
locity to jet velocity. The analysis is compared with experimental data, and the velocity 
profiles at different axial positions are shown to be similar. Centerline velocities follow 
the predicted behavior reasonably well except for the smallest experimental velocity 
ratio; half-radii, an indication of the amount of jet spreading, vary as expected after a 
certain axial distance. The shear-stress distributions are not similar, even at the larg- 
est axial positions for which data are available. In addition, the similarity solution is 
extended to include constant-density mass transfer or heat transfer. 


INTRODUCTION 

The similarity solution for mixing of a turbulent free jet with a quiescent ambient 
stream is well known and is discussed in reference 1. For a circular jet, the flow be- 
comes similar beyond an axial position of about 14 jet radii. Velocity profiles are repre- 
sented by a Gaussian probability function, the centerline velocity decay is inversely pro- 
portional to axial distance, and the half-radius, that is, the position at which the local 
velocity is the average of centerline and coaxial- stream velocities, is proportional to 
axial distance. When the coaxial stream is moving, it is advantageous to use the differ- 
ence between the local axial velocity and the coaxial- stream velocity, rather than the 
local axial velocity alone, as the dependent variable. Radial profiles at a given axial 
position are then expressed as the ratio of this local difference to the maximum (i. e. , 
centerline) difference at that axial position. It was shown by Forstall and Shapiro (ref. 2) 
that when a coaxial stream is moving more slowly than the jet, the flow becomes similar 
after about 8 to 20 radii, depending on the initial ratio of coaxial- stream velocity to jet 



velocity. Velocity profiles can be represented by a cosine, 3/2 power, or error curve. 
Centerline velocity difference was proportional to axial distance, and half-radius was a 
function of axial position and the initial ratio of jet velocity to coaxial- stream velocity. 
The wake behind a disk was investigated by Cooper and Lutzky (ref. 3), who found veloc- 
ity profiles similar beyond 40 radii, centerline velocity difference proportional to the 
-2/3 power of axial position, and half- radius proportional to the 1/3 power of axial posi- 
tion. 

When the coaxial stream is moving faster than the jet, the jet mixing solutions men- 
tioned previously are not applicable. However, in view of the success of these solutions 
in describing velocity profiles, it seems reasonable to ask whether a similarity solution 
exists for turbulent mixing of a jet with a faster moving coaxial stream. If such a solu- 
tion exists, does the real flow exhibit similar velocity profiles, at least in a region of the 
flow field? If so, it only remains to compare the analysis and data to determine the con- 
stants of proportionality that arise in the similarity solution. A similarity solution is ob- 
tained for the turbulent mixing of a jet with a faster moving coaxial stream when the ini- 
tial ratio of coaxial stream velocity to jet velocity is large. The results of the analysis 
are compared with available experimental data that illustrate similarity behavior in part. 


SYMBOLS 


Cp Cg, Cg proportionality constants 

f dimensionless velocity, U/U^ 

g dimensionless mass fraction, Y/Y^ 

o 

h dimensionless shear stress, cr/pU^ 

i dimensionless mass flux, r/pU^Y^ 


K, K j , K 2 , Kg , K 4 constants 

k proportionality constant in eddy viscosity expression 

r radial position 

rj diameter of jet discharge tube 

Sc^ turbulent Schmidt number 

U coaxial- stream velocity minus local axial velocity 

V radial velocity 

W axial velocity 
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x axial position 

Y coaxial- stream mass fraction minus local mass fraction 

e eddy viscosity, kr^yg 1 ^ 

r ? similarity variable, (r/rj )cp 

£ dimensionless axial position, x/rj 

p density 

a turbulent shear stress, pe(9W/9r) 

r turbulent mass flux, (pe/Scp(9Y/9r) 

<p function of axial position 

X dimensionless centerline mass fraction difference, Y^/Y^ 
4? dimensionless centerline velocity difference, U^/Uj 
Subscripts: 

<t centerline 

e coaxial stream 

j initial jet 

1/2 position at which U/U^ =1/2 


ANALYSIS 

Differential Equations and Boundary Conditions 

The equations that will be used to describe turbulent jet mixing are the time- 
averaged continuity equation and the boundary- layer form of the Navier-Stokes equation. 
Molecular transport will be negligible compared with turbulent transport and can be ig- 
nored. The constant-density axisymmetric forms of the equation are (ref. 1) 

-(W)+ii (rV) = 0 (1) 

9x r dr 


and 


pW — (W) + pV — (W) = — — (r or) 
dx dr r dr 


( 2 ) 
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The initial and boundary conditions for the problem are 


W=Wj 0 < r < r j x = 0 


W = W r > r. x = 0 
“ J 


dW 

dr 


= 0, V = 0 r = 0 x > 0 


W = W g r — x > 0 


(3) 

(4) 

(5) 

( 6 ) 


A schematic drawing showing the velocities is given in figure 1. 

When the coaxial- stream velocity is large compared with the initial jet velocity, the 
flow will be more wakelike than jetlike. It is shown in reference 3 that, for wake flows, 
the radial velocity is small far downstream. The approximation that the radial velocity 
is small is used in the present work to simplify the momentum equation; the results will 
thus be restricted to regions far enough downstream that the approximation is valid. It 
will be more convenient to work in terms of a velocity difference 


U = w e - W 


(7) 


than in terms of a velocity. The momentum equation can be made dimensionless with the 
initial velocity difference and jet radius and, when V is small, written as 


( 8 ) 




The boundary conditions are then 



— = 0 


U = 0 




(9) 


JL^oo 2L> 0 

r j r j 


( 10 ) 
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Similarity Solution 


The flow is expected to be similar after a certain distance. In this region, put 

-2-ffo) 

U t 

— = h(r?) 

pvl 

Ux 

U; 


where 


V = — <p(£) 

r j 

and 



r i 


Roman letters are used for functions of the similarity variable r), and Greek letters 
used for functions of the axial variable The momentum equation then becomes 

\ u j / \ d £ d £ d V v d ^? 

with boundary conditions 


— (0) = 0 
dr] 

f(°°) = 0 


( 11 ) 

( 12 ) 

( 13 ) 

( 14 ) 

( 15 ) 
are 

( 16 ) 

( 17 ) 

( 18 ) 
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For downstream, when the coaxial- stream velocity is large compared with the local 
velocity difference (i. e. , W„ » U), equation (8) can be integrated to yield the integral 
momentum equation; if lim (r/r, cr/plL ) = 0, then 


or 



where K is a constant. Equation (20) can be written 
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( 21 ) 


Since the integral is independent of axial position, it is necessary that 





where Kj is a constant. It follows that 


^=2KK 

d? 1 d£ 


so that the momentum equation, when W g » U, becomes 


^ d£ / 2f + v df \ = _kk 1 A. (rjh) 
U- d£ \ dT7/ rj d r\ 


( 22 ) 


(23) 


(24) 


In order for there to be similarity, it is necessary that the terms containing % be 
proportional to the terms containing rj in equation (24). Thus, 
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w 


— ^ = -KICK 9 <p 4 

U j 


( 25 ) 


and 


2 f + 77 — = -i- i — (rjh) 
d 77 K 2 77 dr) 


(26) 


where is a constant. A solution to equation (25) is 


1 U i 

— = 3KK-K,, —L 

3 1 d W 

<p e 


(27) 


The constant of integration has been dropped since it is small compared with the term re- 
tained far downstream. Then, from equation (22), 


\f/ = KKj^ = 


3Kr 


{K 1/2 k 1/2 W e 



-2/3 


(28) 


Equation (26) can be written as 


d ( 77 2 f) = -A- A (rjh) 


d?) 


Kp d?7 


(29) 


and integrated using the boundary condition at 77 = 0 to give 


r/f = — — h 

K 2 


(30) 


The turbulent shear stress can be expressed as 


a = pe 


9W 

dr 


-Piu ^ I 

^ y*u u * 

w 


(31) 


7 



so that 


h 




Thus, equation (30) becomes 



K 2 r j U <t 


€<p 


T]f = 


0 


(32) 


(33) 


It is usually assumed in turbulent jet mixing analyses that eddy viscosity is independ- 
ent of radial position. If Prandtl's postulate that 


e = kr l/2 U t 

is used, the coefficient in equation (33) becomes 

K 2 r j U t _ K 2 1 
ecp k t ? 1/2 


(34) 


(35) 


where ^1/2 = ^ r l/2 //r j^' Since ^i/2> ^at is, the value of 77 at which f = 0. 5, is con- 
stant, equation (33) can be integrated to yield 


f = 



K 2 r j 


2kr 


1 / 2 ^ 



(36) 


Evaluating equation (36) at the half- radius, where f = 0. 5, gives 


1 ^2 r i 2 

In — = -0.693 = £-3- ryf yo 

2 2kr l/2^ 


Thus, the velocity profiles are given by 


f = exp 

- 0. 693 [ v V 

- exp 

- 0.693/^ r Y 


v 1/2/ 

— 


"'csr 

rH 


(37) 


( 38 ) 
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Equation (37) yields 


^ 1/2 


2(0. 693)k 
K 0 


or 


Using equation (22) yields 


: l/2 _ 2(0.693) 1_ 
r j K 2 V 


r iZ? = 2<o.a93 J 1 * 11 ' 172 


r j 


k: 


2 


Equation (22) can be integrated with the help of equation (38) to give 


1 2 _ K cp 

2(0.693) 7?1 / 2 


* 


(39) 


(40) 


(41) 


(42) 


and, since = r i/2^ r j ( f> 


' 1/2 _ 


r j 


2(0. 693)K — 


1/2 


Equations (41) and (43) can be used to eliminate and Kg from equation (28) to give 

(44) 


tt \ “2/3 

if/ = [2(0. 693)]" V 3 /-^ L £ 

(K 1/2 W < 


Then, equation (43) can be written 


= [2(0. 693)1 2//3 (3Kk i 
V W e 


1/3 


(45) 


r i 


I I II 
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Equations (30), (38), and (39) yield the shear-stress distribution 


h = 2(0. 693)k — — 
r l/2 




( 46 ) 


The limiting case when the initial ratio of coaxial- stream velocity to jet-stream ve- 
locity becomes infinite corresponds to wake flow. It is to be noted that the functional 
forms of the expressions for centerline velocity and half-radius (given by eqs. (44) and 
(45)) reduce to those found by Cooper and Lutzky (ref. 3) in their wake-flow investigation. 


Mass Transfer or Heat Transfer 

The analogy between the turbulent transport processes allows the analysis to be ex- 
tended to account for mass transfer or heat transfer as long as density remains constant. 
The same procedure is employed to solve the diffusion equation or energy equation as was 
used to solve the momentum equation; in addition, the solution of the momentum equation 
is used. The dimensionless diffusion equation can be written in terms of a mass fraction 
difference analogous to the velocity difference as 



with boundary conditions 



In the region where the flow is similar put 

— = gfa) (50) 

Y* 
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— - — = H 7 )) (51) 

pU t Y t 

^ = X(9 (52) 

Y ) 

The definitions of 77 and £ are the same as for the solution of the momentum equation, 
namely, 


V = — <P(£) 


r 3 


(53) 


« = 


r i 


The diffusion equation can then be written as 


fe-SpVg*S + X = 

yu. J \ d£ <p d£ dri J ri dri 


(54) 


(55) 


with boundary conditions 

(o) = 0 

dr7 

g(°°) = 0 

Far downstream, when W » U, equation (47) can be integrated to yield 


/ 


Y 2L d JL 1 = K<3 
Y i r i \v 


(56) 

(57) 


(58) 


where Kg is a constant. This equation can be written as 
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(59) 



sn d v = k 3 


or 



(60) 


where is a constant. Then 


^= 2 K„K d cp^ (61) 

d£ a * d£ 


and, when equation (22) is used, equation (55) becomes 


— ^ (2 g + V^) = KKj<p 4 A A (17 i) ( 62 ) 

Uj d| y d r]J 77 dri 

In order for there to be similarity, it is again necessary that the terms containing £ 
be proportional to the terms containing 77 . Thus, in light of equation (25), 

— ^ = -KK-.K^ 4 ( 63 ) 

Uj 


which has already been solved, and 


2g + 77 — = - — — — (r;i) (64) 

d 77 Kg 77 dr? 


This equation can be written as 


— 0? 2 g) = — (f?i) (65) 

d7j Kg dT7 


Integrating and using the boundary condition at 77 = 0 give 
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The turbulent mass flux can be expressed as 



where e is given by equation (34) . Then 

i = e JL dg 

Sc t r j U t d?7 


so that equations (66) and (68) give 


dg K 2 Sc t r j U <t 
d ?7 e<p 


0 


This equation can be integrated to yield 


g = exp 


K 2 Sc t r l U t v 2] 

2ecp 


Using equation (39) gives 


g = exp 


- 0.693 Sc t 



and thus 


g = 



The decay of centerline mass fraction can be obtained when equation (60) is inte- 
grated and equation (71) is used. Thus, using also equation (22), 



X ~ Sc t ^ 


(73) 


For heat transfer, mass-fraction differences are replaced by temperature differences 
and the turbulent Schmidt number is replaced by the turbulent Prandtl number. 


RESULTS AND DISCUSSION 

The functional forms of the variables that describe the turbulent mixing of a jet with 
a faster moving coaxial stream were determined by the analysis. Experimental data are 
required to evaluate the proportionality constants that arise; these constants were ob- 
tained from an experiment (ref. 4) in which time-averaged velocities and Reynolds 
stresses were measured with a constant -temperature hot-film anemometer. The duct 
was large compared with the jet discharge tube in order to minimize the effects of duct 
walls on the flow. Air was used for both the jet and the coaxial stream. Initial ratios of 
coaxial- stream velocity to jet velocity varied from 3. 4 to 39. 5 at a constant Reynolds 

4 

number (based on coaxial- stream velocity and jet discharge tube diameter) of 2x10 . 
Most of the measurements were made within the first 30 jet radii downstream of the jet 
discharge tube. 

Much of this jet mixing data is not far enough downstream to satisfy the restriction 
imposed in the derivation of the similarity solution, namely, that the coaxial- stream ve- 
locity is much larger than the difference between the coaxial- stream velocity and the 
local velocity. However, inspection of the data in figures 2 to 6 shows that the behavior 
predicted by the similarity solution is exhibited at smaller axial positions than would be 
expected. Although the data at the smallest axial positions do no follow this behavior, 
they are included to illustrate that the similarity solution is not everywhere valid. Also, 
the data at the lowest initial ratio of coaxial- stream velocity to jet velocity are included 
to show that the similarity solution is not valid at this low a ratio except possibly at the 
farthest axial position. 

Equations (44) and (45) relate centerline velocity difference and half-radius to rela- 
tive axial position. The jet mixing data were used to determine the constants in these 
equations. The data are plotted in figures 2 and 3, and a line with the slope indicated by 
equation (44) or (45) is drawn through the downstream points. These equations then be- 
come 


^ = 1.88 




(74) 
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The wake-flow data in the downstream region reported (ref. 5) showed that a more rapid 
decay of centerline velocity difference and, correspondingly, a more rapid spreading of 
the mixing region had occurred. The difference in behavior may be due to the boundary 
layer that was present on the jet discharge tube. As a result, the mixing was delayed so 
that a more narrow mixing region was formed. This is consistent with the observation 
that the wake centerline velocity difference data, after the mean eddy pocket, decay as 
the inverse square of axial position before exhibiting the inverse 2/3 power behavior of 
the similarity solution; however, the centerline velocity difference data from the jet mix- 
ing experiments follow the similarity behavior throughout the entire axial region investi- 
gated. 

The relation given by equation (43) shows that half-radius varies inversely as the 
square root of centerline velocity difference and is explicitly independent of axial position. 
Half-radius is plotted as a function of centerline velocity difference in figure 4 and a line 
with a slope of -1/2 is drawn through the downstream points using the constant K deter- 
mined from figures 2 and 3. Equation (43) becomes 


X UA = 0. 60 \f/ 1//2 (76) 

r i 

The downstream data points cluster around the line so determined. The constant is 
nearly the same as that found in reference 5 for wake flow behind a disk. With k deter- 
mined, the eddy viscosity expression can be evaluated and equation (34) becomes 

e = 0.0056 r^ 2 U t (77) 

For comparison with a different physical situation, the proportionality constant for a jet 
mixing with a quiescent ambient stream is about 0. 02 (ref. 1). 

The experimental velocity data for five initial velocity ratios are plotted separately 
in figures 5(a) to (e) along with the line calculated from equation (38); all the data and 
the line calculated from equation (38) are shown in figure 5(f). The agreement at all ve- 
locity ratios is good except at the smallest axial positions. The shear-stress profiles 
are shown in figure 6. The line calculated from equation (46) by using the value of k de- 
termined from figures 2 and 3 is also shown. The profiles are not similar at the largest 
axial distances and fall below the prediction. Predicted shear-stress distributions for 
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turbulent jet mixing with a quiescent ambient stream also lie above experimental data 
(ref. 1). 

Data are not available to determine the proportionality constants for the variations of 
centerline mass fraction or temperature. A reasonable first approximation would be to 
take the turbulent Schmidt or Prandtl number as 0. 7 and use the same proportionality 
constant as centerline velocity difference variation, that is, 1. 88. 


CONCLUSIONS 

A similarity solution for turbulent mixing of a jet with a faster moving coaxial stream 
was derived for large initial ratios of coaxial- stream velocity to jet velocity. The solu- 
tion was compared with data from an experiment in which both the jet and the coaxial 
stream were air. For initial ratios of coaxial- stream velocity to jet velocity from 3. 4 
to 39. 5 and axial distances from 2. 8 to 30 jet radii, the conclusions are as follows: 

1. The radial velocity profiles are similar and follow a Gaussian probability function. 
In terms of a velocity difference U = W g - W, 


where W g is the coaxial- stream velocity, W is the local axial velocity, is the cen- 
terline velocity difference, r is the radial position, and r^yg * s the radial position at 
which U/U^ = 1/2. This even appears to be true upstream of the position where the ex- 
perimental half-radii vary as the similarity solution predicts. 

2. The similarity solution indicates that the centerline velocity can be represented 
as 



Zi-c 

U j \ W e 


where Uj is the initial jet velocity difference; Cj is the proportionality constant, equal 
to 1. 88; x is the axial position; and r^ is the diameter of the jet discharge tube. The 
data follow this behavior reasonably well for the range of axial distances investigated ex- 
cept for the lowest velocity ratio. 

3. The half-radius variation predicted by the similarity solution is 
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1/3 




where Cg is the proportionality constant, equal to 0. 44. The experimental values of 
half-radius cluster around 0. 9 jet radius before following the similarity behavior. 

4. The shear- stress profile predicted by the similarity solution is 


CT = Co — — exp 


pv\ 


1/2 


0.693 



where a is the turbulent shear stress; p is the density; and Cg is the proportionality 
constant, found from the centerline velocity and half-radius variations and equal to 0. 078 
The experimental shear-stress profiles were not similar, even at the largest axial dis- 
tances for which data are available. 

5. The value of the proportionality constant k in Prandtl's hypothesis for eddy vis- 
cosity, e = kr 1 / 2 was determined to be 0. 056. This is about 2^ times the value for 
jet mixing with a quiescent ambient stream. 

6. The analysis was extended to account for constant -density mass transfer or heat 
transfer between a jet and a coaxial stream. The mass-fraction profiles are 



where Y and are the local and centerline values of the coaxial- stream mass frac- 
tion minus the local mass fraction, and Sc^. is the turbulent Schmidt number. Centerline 
mass-fraction variation is given by 



Sc. 


U, 


u ; 


where Yj is the initial jet mass-fraction difference. For heat transfer, mass-fraction 
differences are replaced by temperature differences, and the turbulent Schmidt number is 
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replaced by the turbulent Prandtl number. The experimental data needed to determine 
the constant of proportionality are not available. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 9, 1967, 

122-28-02-16-22. 
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Figure 3. - Axial variation of half-radius. 


initial velocity 
ratio, 
W e fWj 

O 39.5 

□ 28.5 

A 16.0 

V 8.0 

O 3.4 


Eq. (76) 
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Figure 5. - Velocity profiles. 
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(e) Initial velocity ratio, 3. 4. 
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